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Chaotic eigenstates of quantum systems are known to localize on either side of a classical partial transport bar¬ 
rier if the flux connecting the two sides is quantum mechanically not resolved due to Heisenberg’s uncertainty. 
Surprisingly, in open systems with escape chaotic resonance states can localize even if the flux is quantum 
mechanically resolved. We explain this using the concept of conditionally invariant measures from classical 
dynamical systems by introducing a new quantum mechanically relevant class of such fractal measures. We 
numerically find quantum-to-classical correspondence for localization transitions depending on the openness of 
the system and on the decay rate of resonance states. 
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Localization of quantum eigenstates and wave packets is of 
fundamental importance for the physics of transport and ap¬ 
pears for a variety of reasons, e.g., strong localization due to 
disorder HI], weak localization due to time-reversal symme¬ 
try [ 2 ], localized edge states due to topological protection [3j], 
or localization due to classically restrictive phase-space struc¬ 
tures 0 • In the latter case, the localization can originate from 
impenetrable barriers of regular motion or partial transport 
barriers with a small transmission given by a flux within 
a chaotic region [4-Ll]. Such partial barriers are ubiquitous 
in the chaotic region of generic two degree-of-freedom Hamil¬ 


tonian systems [5!, [ 6 ], [90 and a universal localization transition 
was found id. Chaotic eigenstates of the system typically 
localize on either side of a partial barrier if the transmission 
region is quantum mechanically not resolved, i.e., if the clas¬ 
sical flux <f> across the partial barrier is much smaller than the 
size h of Planck’s cell (T> <C h). If the transmission region 
is quantum mechanically resolved (h <F), eigenstates are 
equipartitioned in the chaotic component, thereby ignoring the 
presence of the partial barrier. 

In contrast, in open Hamiltonian systems which allow for 
escape Il3l423l. chaotic resonance states exhibit localization 
in the presence of a partial barrier surprisingly even in the 
semiclassical regime (h <C <F) 10]. Such a localized state is 
shown in Fig. [U upper right, by its Husimi phase-space rep¬ 
resentation. This demonstrates that in open systems the influ¬ 
ence of partial barriers on localization and transport proper¬ 
ties is even more substantial than in closed systems. A thor¬ 
ough understanding of this localization phenomenon remains 
open, so far. A prominent application are optical microcavi¬ 
ties, where the emission patterns are governed by the localiza¬ 
tion of eigenmodes 1251-031. For their design, it is particularly 
important to know whether a partial barrier is desired to en¬ 
hance localization or whether it should be avoided. The local¬ 
ization phenomenon may also have relevance in many other 
areas of physics, such as transport through quantum dots 0 , 
ionization of driven Rydberg atoms [ 35 ], and microwave cav¬ 
ities [36]. 

Since the localization appears in a semiclassical regime 
(h T>), one may wonder if it has a classical origin. Thus one 



FIG. 1. (color online) Weight ||Pi / 0 7 || 2 (symbols) of resonance 
states in region A± vs ratio of size \Q\ of opening and flux <h across a 
partial barrier for different parameters of the partial-barrier standard 
map (16 < <k//i, |fi| /h < 2048; |Ai| = 0.5; h = 1/6000). Weight 
of state with 7 closest to 7 na t (red points) and averaged over states 
with decay rates 7 G [ 7 nat/1 -1,1.1 7 nat] (black crosses). This is com¬ 
pared to the natural CIM p na t (Ai) [Eq. ([4]), solid green line]. Inset: 
Phase space of the partial-barrier map, illustrating regions A±, A 2 on 
either side of the partial barrier (solid magenta line) with exchanging 
regions <3>i, $ 2 , and opening Q. Upper panels: Husimi represen¬ 
tation of typical resonance states with 7 « 7 na t for h = 1 / 1000 , 
&/h = 20 , and two values |U|/<f> indicated by arrows. 


needs the classical counterpart of a quantum resonance state. 
This is given in the field of open dynamical systems HI 53- 
44| by a conditionally invariant measure (CIM). It is invariant 
under time evolution up to an exponential decay with rate 7 . 
The asymptotic decay of generic initial phase-space distribu¬ 
tions leads to the so-called natural CIM /i nat with decay rate 
7nat- The quantum-mechanical relevance of /i na t is shown in 
I13- H51 I2A Ull. Note that the steady probability distribution 
introduced in the context of optical microcavities [ 0 ] corre¬ 
sponds to /inat- The natural CIM /x na t for the single decay rate 
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FIG. 2. (color online) Weight ||Pip 7 || 2 (red points) of resonance 
states in region A\ vs decay rate 7 for the partial-barrier standard 
map (&/h = 64; \Q\/h = 1024; |Ai| = 0.5; h = 1/6000). This 
is compared to the 7 -natural CIM p^(Ai) [Eq. 0), solid green line]. 
Upper panels: Husimi representation of typical long-lived (left) and 
short-lived (right) resonance state for h — 1/1000 with 7 values 
indicated by arrows. 


7 nat, however, cannot be the classical counterpart for all quan¬ 
tum resonance states as they have a wide range of decay rates 
(see, e.g., Fig. [2]). Exceptional CIMs with decay rate 7 differ¬ 
ent from 7nat have been discussed Hi. In fact, for each 7 
one can construct infinitely many CIMs. It is an open question 
which of these CIMs correspond to quantum resonance states 
for arbitrary 7. To answer this question one has to go beyond 
the important results of Ref. 118] which relate the total weight 
of a resonance state on each forward escaping set to its decay 
rate. 


In this Letter, we introduce the quantum mechanically rel¬ 
evant class of CIMs. Their localization explains the localiza¬ 
tion of chaotic resonance states in the presence of a partial 
barrier. In particular, we find (i) a transition from equiparti- 
tion to localization when opening the system, Fig. [TJ and (ii) 
a transition from localization on one side of the partial bar¬ 
rier to localization on the other side for resonance states with 
increasing decay rate, Fig. [2 We numerically demonstrate 
quantum-to-classical correspondence for a designed partial- 
barrier map and the generic standard map. 


Partial-barrier map .— We design a chaotic model map with 
a single partial barrier (similar to Ref. (24]), which allows for 
numerically varying the flux across the partial barrier and for 
deriving the classical localization, Eq. The partial-barrier 
map T = MoFoOisa composition of three maps: the map 
M describes the unconnected chaotic dynamics within two 
regions, A They decompose the phase space T = [0,1) x 
[0,1) into A\ = [0, | Ai|) x [0,1) and its complement A 2 = 
r \ Ai\ see the inset in Fig. [2 where \Ai\ denotes the area 
of A 1 and by normalization, |A 2 | = 1 — \Ai\. The map E 


induces a flux T> between A\ and A 2 by exchanging regions 
<&k C Ak with \<$>k\ = The map O opens the system by the 
absorbing region $2, which is contained in region A\. 

We introduce two different dynamics for M. For the nu¬ 
merical analysis, we use the generic standard map 0 on 
the torus in symmetrized form, q t + 1 = q t + p*, pt +1 = 
Pt + v (<h+ 1 ) with p* = p t + v(q t ) for v(q) = £ sin(27rg) 
acting individually on each of the regions Ak after appropriate 
rescaling. We fix n = 10 where the standard map displays a 
fully chaotic phase space. For analytical considerations, we 
use the ternary Baker map in each region Ak , as illustrated in 
Fig. 0a), which allows for the derivation of Eq. @. We refer 
to the corresponding maps T as partial-barrier standard map 
and partial-barrier Baker map , respectively. 

Quantum localization transitions. —Let us consider the 
quantization U of the partial-barrier standard map T. From 
the eigenvalue problem for U , 

u^ (i) 

we numerically compute the decay rates 7, describing the tem¬ 
poral decay of the norm, ||LPt /> 7 || 2 = e _7t , and the corre¬ 
sponding resonance states t /) 7 (the phase 6 is not relevant in 
the following). The absolute weight of t /? 7 in region A\ is 
given by ||Pit/) 7 || 2 , where Pi denotes the projection onto the 
subspace associated to A\. We observe (i) a transition from 
equipartition to localization on A 2 for increasing size | O | of 
the opening, see Fig. [T] and (ii) a transition from localization 
on A 2 to localization on Ai for increasing 7, see Fig. [21 Tran¬ 
sition (i) is surprising as localization occurs for h <C T>, where 
in the closed system all eigenstates are equipartitioned ill. 
Transition (ii) shows that in open systems the localization de¬ 
pends on the decay rate 7. 

In Fig. [TJ we focus on resonances with decay rate 7 ~ 7 nat , 
which describes the decay of typical long-lived resonance 
states in the semiclassical limit. We find transition (i) from 
equipartition, ||Pit /) 7 || 2 = |Ai|, for \Q\ T> to localization 
on A 2 for \Q\ <4> for various values of 3>//i and \Q\/h. 

The transition is universal with the scaling parameter |fi|/<£. 
Moreover, this even holds for individual states without aver¬ 
aging (red dots). We stress that this localization transition in 
the open system occurs even though &/h > 10 , where in the 
closed system all eigenstates are equipartitioned [ 12 ]. 

In Fig. [2 we fix the parameters such that \Q\ T>, for 

which the long-lived resonance states localize on A 2 , and 
show the 7 dependence of the weights ||Pit /) 7 || 2 for all res¬ 
onance states. We find transition (ii) from resonance states 
which localize on A 2 for small 7 to resonance states which lo¬ 
calize on A\ for large 7, including equipartitioned resonance 
states in between. 

The fact that both transitions (i) and (ii) occur for h <f> 
suggests that the localization transitions could be of classical 
origin. Furthermore, from the point of view of decaying clas¬ 
sical distributions the observed transitions qualitatively seem 
to be rather intuitive: in Fig. [2 for a larger size of the opening 
one has less weight in region A 1 . In Fig. [2 a larger weight 
in Ai corresponds to a larger decay rate. For a quantitative 
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FIG. 3. (color online) (a) Illustration of the partial-barrier Baker 
map T = M o E o O. Magenta line indicates partial barrier and 
gray shaded region marks the opening (left and central) and image of 
opening (right), (b) Backward trapped set (dark horizontal stripes) 
and forward escaping sets Cl (gray), T - 1 (Q) (yellow), T~ 2 (Cl) (or¬ 
ange), and T~ 3 (Cl) (red), (c) Natural CIM integrated over boxes of 
size 3 -3 in the p direction, (d), (e) Approximation of 7 -natural CIMs 
by truncation of Eq. © to n < 2 for 7 / 7 na t. 


description, however, one needs to find the quantum mechan¬ 
ically relevant class of CIMs. 

Classical localization. —A conditionally invariant measure 
(CIM) /i 7 is defined by 

/x 7 (T- 1 (X))=e-> 7 (X), ( 2 ) 

for each measurable subset X of phase space. It is invariant 
under the classical iterative dynamics T of the open system 
up to an exponential decay with rate 7 . Equation ([2]) states 
that the measure /i 7 (T _ 1 (X)) of the set T~ 1 (X) that will 
be mapped to X is smaller than p 7 (X) by the factor e -7 . 
These measures must be zero on the iterates of the opening Cl. 
Thus, the support of /i 7 is the fractal backward trapped set 
Tb [horizontal black stripes in Fig. Ob)], that is the set of 
points in phase space which do not escape under backward 
time evolution. Particularly important is the natural CIM /i nat , 
see Fig. 0c), which is constant on its support [because of in¬ 
tegration over boxes in Fig. 0c) one finds two nonzero box 
measures]. 

We now generalize fi nat to a CIM /i 7 of arbitrary decay rate 
7 , which we call 7 - natural CIM. To this end, we use a con¬ 
struction of CIMs [40, [4 it] where one starts with an arbitrary 
probability measure on the intersection Cl D Tb of the open¬ 
ing Cl with the backward trapped set Tb. By propagating this 
measure backwards to all forward escaping sets T~ n (Cl ) [ver¬ 
tical colored stripes in Fig. 0b)] and appropriate scaling [re¬ 


specting the decay rate 7 , Eq. ©] one obtains a CIM. Here, 
we choose the simplest measure on Cl n Tb, given by /i na t- 
This choice of a measure, which is constant on its support, is 
quantum mechanically motivated in analogy to quantum er- 
godicity for closed fully chaotic systems, where eigenstates in 
the semiclassical limit approach the constant invariant mea¬ 
sure [£ 6 , 47]. This choice leads to the 7 -natural CIM 

00 

fh(X) e (7 “'- 7) >na,(X n T- n (fi)), (3) 

n =0 


with normalization J\f = (1 — e _7 )/(l — e _7nat ). This se¬ 
ries multiplies /i nat in each forward escaping set T~ n (Cl) by 
an appropriate factor which imposes the overall decay rate 7 
according to Eq. (0). Two examples of 7 -natural CIMs for the 
partial-barrier Baker map are shown in Figs. [3jd) and 0e). 
The measure is constant on D Tb for each n G No. 

With increasing n, this constant is decreasing (increasing) for 
7 > 7 nat (7 < 7 nat); in particular, short-lived measures /i 7 
have more weight in the opening. Note that the idea underly¬ 
ing Eq. ([3]) was used without the notion of CIMs in Ref. [18] 
for sets X = T~ n (Cl ) for systems without a partial barrier. 
Moreover, note that the 7 -natural CIMs are solutions of the 
exact Perron-Frobenius operator (which is not available), but 
cannot be obtained from finite-dimensional approximations. 
Therefore, they have to be constructed directly in phase space. 

We find as our main result on the classical localization of 
/i 7 due to a partial barrier that the weight of /i 7 on each side 
of the partial barrier is given by 0 


/i 7 (Ai| — 


1 — Cry 


(4) 


and // 7 (A. 2 ) = 1 — /x 7 (Al), with 

c 7 = (l — e 7_7nat ) (l — e _7nat ) 


j^ij j^gj 


(5) 


The values for /i n at(^i) and 7 nat follow from the longest-lived 
eigenstate of the eigenvalue problem 


f /^nat(-4-l) N \ 
\/fii a t(-4.2) / 


— g 7nat 


f Mnat( 7 4l) N \ 
\/fii a t(-4-2) ) 


( 6 ) 


where F nat denotes the transition matrix between A\ and A 2 
for the one-step propagation of p nat (see Ref. (49] for approx¬ 
imations of the Perron-Frobenius operator). In general, F nat 
may be obtained numerically or it may be approximated by 
assuming a uniform distribution for /i nat , 


/l-dni + ^/IAxi <71^21 \ 

V */\M i-$/\a 2 \J- 


(7) 


This turns out to be quite a good approximation even for frac¬ 
tal /inat and it is exact for the partial-barrier Baker map. 

Quantum-to-classical correspondence. —Figure |T] (green 
line) shows the classical localization p 7 (Ai), Eq. (0]), for 
7 = 7nat C7 = 0 and /i 7 = /i na t)> using the approxi¬ 
mation Eq. ©. When increasing the size |0| of the opening, 
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FIG. 4. (color online) Weight ||Pi ^ 7 || 2 (red points) of resonance 
states in region 4i vs decay rate 7 for standard map at k = 
2.9, with |Ai| « 0.6664, |A 2 | « 0.2061, <f> « 0.0126, |Q| = 
0.1, and h = 1/10000. This is compared to the 7 -natural CIM 
/i 7 (Ai), either by direct numerical computation [48] (solid green 
line), using Eq. © and computing F nat numerically (dashed green 
line), or by using approximation Eq. (dotted green line). Inset: 
Phase space of the standard map with regular and chaotic regions, 
illustrating regions A 1 (medium gray shaded), A 2 (light gray shaded) 
on either side of the main partial barrier (thick solid magenta line) 
with exchanging regions 4>i, <f> 2 , and opening Q (dark gray shaded). 
Upper panels: Husimi representation of typical long-lived (left) and 
short-lived (right) resonance states for h = 1/1000 with 7 values 
indicated by arrows. 


we find a transition from equipartition for |Q| 4> to local¬ 
ization for \Q\ T>. The only scaling parameters are |Q|/4> 

and |Ai |/|A 2 |. We find very good agreement of the classical 
localization measure with the localization of the quantum res¬ 
onance states. Note that for 7 7 ^ 7 nat , the localization depends 
on all parameters |Q|, <E, and | A\ |. 

Figure [2] (green line) shows Eq. 0 as a function of 7 . 
The classical localization measure /i 7 (Ai) monotonically in¬ 
creases with 7 ; i.e., the faster the decay, the larger is the 
weight in region A\ with the opening Q. In the limit 7 oc 
one finds /x 7 (Ai) = 1 , and in fact, all the weight is in the 
opening Q. In the limit 7^0 one finds a small constant 
lio(Ai) > 0 ; i.e., even though most of the weight is in A 2 
there is always a small contribution in A\ due to the exchange 
between A 1 and A 2 . Again, we find very good agreement be¬ 
tween the classical localization measure and the localization 
of the quantum resonance states for all decay rates 7 . Note 
that quantum-to-classical correspondence is also confirmed 
for | Ai| 7 ^ IA 2 I (not shown). 

Do these results for the partial-barrier map generalize to 
generic systems? In Fig. 01 we show for the standard map at 
k = 2.9, where it has a mixed phase space, that the local¬ 
ization of the chaotic resonance states on region A 1 , which 
contains the opening, increases as a function of 7 . Qual¬ 


itatively, we find the same localization behavior as for the 
partial-barrier standard map in Fig. [2 Quantitatively, it is 
well described by the classical localization of /i 7 , which is 
determined numerically 0 . Also the analytical prediction, 
Eq. 0, works reasonably well. Overall, Figs. [U [2 and 0| 
demonstrate quantum-to-classical correspondence for the lo¬ 
calization of chaotic resonance states in open systems due to 
a partial barrier. 

Outlook .—We see the following future challenges: 
(a) While in this work we concentrate on the weights 
on either side of a partial barrier one should verify the 
quantum-to-classical correspondence for the fine-structure 
of chaotic resonance states to 7 -natural CIMs. (b) Which 
deviations arise when approaching the quantum regime 
of ft « |fi|? (c) Is the new class of 7 -natural CIMs, 

which is quantum mechanically motivated, of relevance 
also in classical dynamical systems? (d) Is it possible to 
predict which quantum mechanical decay rates 7 occur in the 
presence of a partial barrier including their distribution, as 


it is known for fully chaotic systems [20, 5y, 51]? (e) The 


present work explains the localization of resonance states 
which have been used to derive the hierarchical fractal Weyl 
laws Q for a hierarchy of partial barriers. Now it is possible 
to discuss whether these laws survive in the semiclassical 
limit, (f) We see direct applications to mode coupling in 
optical microcavities [52] and in recently studied parity-time 
symmetric systems 1531 [54j], where instead of a partial barrier 
one has coupled symmetry-related subspaces. 
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stimulating discussions, and acknowledge financial support 
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SUPPLEMENTAL MATERIAL 

Classical derivation .—In order to derive Eq. we focus 
on the localization properties of /i 7 with respect to the partial 
barrier, restricting ourselves to the partial-barrier Baker map 
in the following. The generalization to other systems will be 
discussed at the end. 

The localization of /i 7 is described by its weights /i 7 (A/~) 
on either side of the partial barrier. In virtue of Eq. (3), 
we only have to study /i n at(^/c H T~ n (fl)) in more detail 
to compute /i 7 (A/~). We find that the natural measure of 
Ak fl T~ n (fl) is proportional to its relative area inside Ak, 

Mnat (A k fl = Mna t(A k ) • l Xl}l (SI) 

\ A k\ 

This follows from the fact that the forward escaping sets 
T~ n (fl) decompose the backward trapped set Tb in the un¬ 
stable (horizontal) direction, on which /i nat is uniformly dis¬ 
tributed within A\ and A2 individually, see Fig. 3(c). 

The distribution of the opening fl over phase space under 
backward time evolution, which enters Eq. (fSTh in terms of 
|Ak fl T~ n (fl)\, follows from 


is indicated by the numerical findings for the partial-barrier 
standard map, Fig. 2, as well as the standard map, Fig. 4, com¬ 
pared to Eq. (4). The main reason for the deviations of Eq. (4) 
for a generic system is that Eq. (ED is not valid in general. 
For its generalization, one has to revise the argument in the 
natural coordinates of the stable and unstable manifolds. 

Standard map .—In the following, we explain how to de¬ 
termine /i 7 for a generic system like the standard map nu¬ 
merically. First, one has to approximate (the chaotic part of) 
the backward trapped set T\>. To this end, one may define a 
uniform grid of N gr id points in phase space of which one has 
to discard points which leave the system within Wi ter itera¬ 
tions of the map T in backward time direction. For Fig. 4, we 
choose TVgnd = 10 6 and TViter = 50. Points within a regular 
phase-space region should be omitted manually. The remain¬ 
ing points provide the finite-time approximation r{J um of Tb 
and need to be classified by their forward escaping times. Fol¬ 
lowing Ref. [18], we associate the /i 7 measure e _7n (l — e -7 ) 
to set T~ n (fl) fl r{J um . Finally, assuming equidistribution for 
the points in T~ n (fl) fl r{J um , we find 

H (X H T~ n (fl)) « f n (X) e _7n (l - e“ 7 ), (S4) 


= F n /W 
\\A 2 nT- n (n)\J nat ^o )' 


(S2) 


where F nat denotes the transition matrix between A\ and A 2 
for the one-step propagation of p mt . Note that the transition 
matrix for the backward time evolution of fl is given by ^nat 
itself. This relation can be interpreted using Fig. 3(b): In the 
beginning, fl (gray vertical stripe) is supported on 4i. In the 
next step, T~ 1 (fl ) (yellow) splits into equal parts on A\ and 
A 2 . Afterwards, T~ 2 (fl) (orange) contributes two stripes to 
A 1 and three to A 2 . 

Inserting the relations (fSTI) and (IS2b in Eq. (3), and using 
Neumann’s series, we obtain 


llry^Afo) - J\f 


\ A k\ 


(a -r-' 1 



(S3) 

Our result on the classical localization, Eq. (4), follows from 
Eq. (IS3l) after some algebra. 

We believe that this derivation for the partial-barrier Baker 
map can be generalized to generic dynamical systems. This 


for each (measurable) subset X of phase space and 


= #(xnr-»(»)nrr m ) 
Jn{ # (T~ n (n) n r“ m ) 


(S5) 


where # denotes the cardinality of a set. Using p^ 7 (X) = 
51 T~ n (fl )) we have a numerical estimate for the 
7 -natural CIM /i 7 . As the sample r{J um is only finite the series 
will terminate and the numerically approximated measure is 
not perfectly normalized. This method is not appropriate for 
exceedingly small 7 since the weight on forward escaping sets 
T~ n (fl) nr{J um with large escape times n increases while they 
are approximated by a few points only. Hence, we restrict the 
method to 7 > 7 na t in Fig. 4. 

Quantum mechanically, owing to the mixed phase space of 
a generic system, we discard all regular and deeper hierarchi¬ 
cal states having less than 50% of their weight within A\ and 
A 2 . As some of the remaining chaotic resonance states still 
have a significant contribution outside of A± U A 2 , we renor¬ 
malize them such that HPiV^II 2 + ||F > 2 ' 0 7 1| 2 = 1. 







